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The effects of the finite size of a liquid drop undergoing a phase transition are described in terms 
of the complement, the largest (but mesoscopic) drop representing the liquid in equilibrium with the 
vapor. Vapor cluster concentrations, pressure and density from fixed mean density lattice gas (Ising) 
calculations are explained in terms of the complement generalization of Fisher's model. Accounting 
for this finite size effect is important for extracting the infinite nuclear matter phase diagram from 
experimental data. 
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PACS numbers: 

Finite size effects are essential in the study of nuclei 
and other mesoscopic systems for opposite, but comple- 
mentary reasons. In modern cluster physics, the prob- 
lem of finite size arises when attempts are made to relate 
known properties of the infinite system to cluster prop- 
erties brought to light by experiment For nuclear 
physics, the problem is the opposite: finite size effects 
dominates the physics at all excitations and the chal- 
lenge is generalize specific properties of a drop (nucleus) 
to a description of uncharged, symmetric infinite nuclear 
matter. This goal has been achieved already for cold nu- 
clei by the liquid drop model. Finite size effects are also 
encountered in nuclear physics in efforts to generate a 
liquid- vapor phase diagram from heat capacity measure- 
ments (2|] and fragment distributions (3]. 

We present a general approach to deal with finite size 
effects in phase transitions and illustrate it for liquid- 
vapor phase coexistence. A dilute, nearly ideal vapor 
phase is in equilibrium with a denser liquid phase; finite- 
ness is realized when liquid phase is a finite drop. A finite 
drop's vapor pressure is typically calculated by including 
the surface energy in the molar vaporization enthalpy 
0, 13 . We introduce here the concept of the complement 
to extend and quantify finite size effects down to drops 
as small as atomic nuclei. We generalize Fisher's model 

, deriving an expression for cluster concentrations of a 
vapor in equilibrium with a finite drop and recover from 
it the Gibbs-Thonison formulae 8] . We demonstrate our 
approach with the lattice gas (Ising) model. 

The complement method consists of evaluating the free 
energy change occurring when a cluster moves from one 
phase to another. For a finite liquid drop in equilibrium 
with its vapor, this is done by virtually transfering a 
cluster from the liquid drop to the vapor and evaluating 
the energy and entropy changes associated with both the 
vapor cluster and the residual liquid drop (complement). 
This method can be generalized to incorporate energy 
terms common in the nuclear case: symmetry. Coulomb 
(with caution j^) and angular momentum energies. 

In the framework of physical cluster theories of non- 



ideal vapors (which assume the monomer-monomer in- 
teraction is exhausted by the formation of clusters), clus- 
ters behave ideally and are independent of each other. 
The complement method is based upon this indepen- 
dence. Physical cluster theories state that the concen- 
trations of vapor clusters of A constituents nA{T) de- 
pend on the free energy of cluster formation AG/i(T) = 
AEa{T)-TI^Sa{T). The epigon of physical cluster the- 
ories is Fisher's model which writes, at coexistence, 
/\E = cqA" and AS'a(T) = ffA'" - rlnA. Thus 



nA{T) = exp 



AGa(T) 



T 



= qaA exp - 



T 



(1) 



where qq is a normalization, r is Fisher's topological ex- 
ponent, Co is the surface energy coefficient, a is the sur- 
face to volume exponent, and e = (Tc — T)/Tc. The 
leading term in A5yi(T), proportional to A"^ , permits 
the vanishing of the cluster free energy at a T = Tc inde- 
pendent of size. Equation |^ (and the extension below) 
is valid only at phase coexistence for T <Tc. The direct 
physical interpretation of the parameters in AG a (T) and 
their application to the nuclear case is the reason for this 
choice here, despite its limitations 9]. 

We generalize eq. valid for infinite liquid- vapor 
equilibrium, to the case of a vapor in equilibrium with 
a finite liquid drop. For each vapor cluster we can per- 
form the gedanken experiment of extracting it from the 
liquid, determining entropy and energy changes of the 
drop and cluster, and then putting it back into the liq- 
uid (the equilibrium condition) , as if, according to physi- 
cal cluster theories, no other clusters existed. Fisher's 
expressions for AEa{T) and ASa{T) can be written 
for a drop of size in equilibrium with its vapor as 
(Ai-Ar-A'^^] and A^^(r) = 
-A'^]-T In [A{Ad - A) /Ad] giving 



AEAiT) = Co [A'' 
f [A'^ + {Ad - AY 



riAiT) 



A {Ad -A) 
Ad 



expi-— K + (Ad-A)^ 



]}• (2) 
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FIG. 1: On the left (right), top to bottom: the d = 2 (3) 
liquid drop size A^, specific heat Cv- See text for details. 



The free energy cost of complement {A^ — A) formation 
is determined just as the free energy cost of cluster (A) 
formation is determined. The resulting expression Eq. (0) 
reduces to Eq. |^ when Ad tends to infinity and contains 
exactly the same parameters. We can rewrite Eq. (jJJ as 



(3) 



with rf^{T) given by Eq. (0). The finite size of the 
drop acts as an effective chemical potential, A/Xfs = 
- {coe [(^d - Af - Al] - Tr In [(Ad - A) /Ad]} /A, in- 
creasing the vapor pressure 

In order to quantitatively demonstrate this method, we 
apply it to the canonical lattice gas (Ising) model [islll^ 
with a fixed number of up spins, i.e. a fixed mean oc- 
cupation density pgxed lattice gas (eqtiivalently, a fixed 
magnetization Mgxed Ising model) [l^- Up spins rep- 
resent particles of the fluid forming monomers, dimers, 
drops etc. Down spins are empty space; the lattice is 
the container enclosing the fluid. We chose large lattices 
with periodic boundary conditions to minimize finite lat- 
tice effects, irrelevant to our study. For d = 2 (3) we 
used a square (simple cubic) lattice of side L = 80 (25) 
which leads to a shift in of < 0.5% [H [13 (< 0.5% 
0,Q|). The Mfixed calculations were performed accord- 
ing to ref. 0. For every (T, pgxed) over 10^ thermal- 
ized realizations were generated to produce the cluster 
concentrations. We performed Mfree calculations for the 
same lattices as a benchmark in order to differentiate 
effects of a finite lattice from those of a finite drop. 




FIG. 2: [Color online] The cluster concentrations of the d = 2, 
1/ = 80 periodic boundary condition square lattice for: AffVcc 
calculations (top); Mflxed calculations with no complement 
(middle); Mflxed calculations with the complement (bottom). 



For the Affixed calculations at T = 0, the up spins ag- 
gregate into a single liquid drop in a vacuum: the ground 
state. At higher temperatures, the vacuum is filled with 
a vapor made of up spin clusters. Clusters in the vapor 
were identified via the Coniglio-Klein algorithm |2oj to 
insure that their behavior is physical (i.e. cluster con- 
centrations return Ising critical exponents and not per- 
colation exponents). The largest cluster represents the 
liquid drop and is identified geometrically (all like spin 
nearest neighbors bonded) in order to capture the skin 
thickness associated with liquid drops |l9|. Our choices 
of ground state liquid drops A'^ (shown in Fig^ insure 
that the ground state is approximately square (cubical) 
for d = 2 (3). Due to periodic boundary conditions the 
ground state shape changes with Affixed 

Figurenshows that as T increases, the drop's size A^ 
decreases from A'^; the evaporating drop fills the con- 
tainer with vapor. At a temperature Tx, corresponding 
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FIG. 3: [Color online] Same as Fig.Hbut for the d = 3, L = 25 
periodic boundary condition simple cubic lattice. 



to the end of two phase coexistence, falls quickly. The 
value of Tx varies with Aj. For T < Tx the specific heat 
Cv (measuring spin-spin interaction energy fluctuations 
only) agrees approximately with the Mfroc results (solid 
lines in Fig. for all until fluctuations in A^ at ^ Tx 
produce a Cy greater than that of the Mfrcc results. As 
T increases further, Cy decreases. 

To evaluate the efflcacy of the complement, we exam- 
ine the scaled cluster concentrations for our calculations: 
nA{T)/qoA~''^ vs. coA'^e/T. For Mftcc calculations it has 
been shown that this scaling collapses concentrations of 
clusters over a wide range in A and T [^. Finite size 
liquid drop effects will be manifested in the cluster con- 
centrations of the Affixed calculations which should scale 
better according to Eq. ^ than to Eq. 

Only clusters of ^ > 9 are included in the Mfrec fits for 
the d = 2 calculations. This is because only large clus- 
ters obey Fisher's ansatz for the cluster surface energy: 
Ea = cqA"^ . Small clusters are dominated by geometrical 



TABLE I: Fit results for M^cc calculations 



- Onsager this work theoretical values this work 

- d=2 d = 2 d = 3 d = 3 
-L-^oo L = 80 L -* oo L = 25 

xl - 4.7 - 1972.2 

Tc 2.26915 2.283 ± 0.004 4.51152 ± 0.00004 4.533 ± 0.002 

Co >8 8.6 ±0.2 > 12 12.63 ±0.04 

a 8/15 0.56 ±0.01 0.63946 ± 0.0008 0.725 ±0.003 

r 31/15 2.071 ±0.002 2.209 ± 0.006 2.255 ± 0.001 



TABLE IL xt results for Mflxcd calculations 



yld {d = 2) 640 


320 


160 


xt w/o complement 10.3 


10.6 


18.2 


xt w complement 1.7 


1.9 


4.8 


A'i {d = 3) 468 


234 


117 


xt w/o complement 14 825.9 


7 938.6 


3 516.0 


xt w complement 1 553.4 


838.6 


258.1 



shell effects [i^- For the c? = 3 Mfixod calculations, large 
clusters are very rare, so clusters oi A>2 are included in 
our analysis. Thus, the magnitudes of the xt values from 
the c? = 3 calculations are due to the clusters analyzed 
not following closely Fisher's ansatz. 

In the thermodynamic limit, the highest temperature 
admitted by Eq. Q or |(5J) is the temperature at which 
the system leaves coexistence. For the Mfree calculations 
this occurs at T = while for the iWfixed calculations 
this occurs at T « Tx- However, due to the small size 
of our drops, fluctuations grow large before Tx and we 
consider only T < OMTx (O.TST^ ) for d = 2 (3). 

To make a comparison between the scaling achieved 
with the Mfrcc clusters and with the Mfixod clusters, we fit 
the Mfrcc clusters with Eq. Q with the free parameters 
Tc, Co, (J and r; go — C{t ~ Results are given in 

Tableland shown in top panels of figures [5] and O The 
values of the Tc and r returned by this procedure are 
within 1% of their established values. The value of co 
shows that all clusters are not perfect squares (cubes) in 
d = 2 (3) for which co = 8 (12) The values of a 

are within 5% (15%) of their established values for d = 2 
(3) . This level of inaccuracy arises from shell effects . 

Next we calculate xt the Mfjxed clusters using 
Eq. (0) and Eq. with parameters fixed to the Tabled 
values for the infinite system. This procedure frees us as 
much as possible from the drawbacks of Fisher's model 
so we can concentrate on the effect of the complement. 
Results are given in Tabic |n] and shown in the middle 
panels (without complement: concentrations scaled as 
nA(T)/qoA~'^) and the bottom panels (with complement: 
concentrations scaled as n^(r)/go^~^ exp [AAufs/T]) of 
figures 121 and 13 The Mffxcd xt values for the calculation 
with the complement are an order of magnitude smaller 



4 



X 
P 

P 

-a 

< 



P 

P 
< 



1.2 
1 

0.8 
0.6 
0.4 
0.2 

1.2 
1 

0.8 
0.6 
0.4 
0.2 



1 ^ ' ' 1 ' ' ' 1 ' ' ' 1 

- AO =160 : 
AO = 320 : 
: - A0 = 640 : 
d = 2: 


.1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 l. 

; ' ao=ii7- 

; AO =234; 

- aO = 468- 
d = 3 


L L i i 1 i i i 1 i i i 1 ; 
: d = 2: 

■ 1 , , , 1 , , , 1 , , , 1 " 


d = 3- 



1.2 



1.4 



1.6 0.75 
T 



1 1.25 1.5 



175 
T 



FIG. 4: Left (right) the normalized pressure and density of a 
vapor in coexistence with a drop Ad for d — 2 (3). Thin [thick] 
lines show no complement results, X = Y — 1 [complement 
results, X and Y from Eqs. ^ and Q]. 



than the results without the complement and the data 
collapse is better. 

The middle and bottom panels of figures El and 01 show 
the major point of this paper: generalizing Fisher's model 
with the complement accounts for the finiteness of the 
liquid. The middle panels of figures [21 and show that 
the Mfixcd calculations do not scale as their Mfrcc coun- 
terparts. The bottom panels of figures jH and O show that 
when the complement effect is taken into account, the 
Mflxcd calculations scale as their Mf^oo counterparts. 

We now confront the integrated quantities, pres- 
sure p{Ad,T) = TJ2a''^a{T) and density p{Ad,T) = 
J2a''^a{T)A . For Ad ^ A, expanding A^fg gives 



A/xfe = 1 + A 

This leads to 

p{Ad,T) « poo(7')exp 
= PooiT)X 

and 

p{Ad,T) « /9oo(T)cxp 
= PooiT)Y. 
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TA^-" 



(4) 
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For a vapor of monomers as A^ ^ r equations (0) and 
© yield the Gibbs-Thomson formulae jsj . 

Figure^shows the behavior of p{Ad,T) and p(yld,T) 
for the Mgxed calculations compared to the bulk results. 
To free ourselves from finite lattice size effects pac,{T), 
Poo (T) and (T) are determined from the Mf^ee calcu- 
lations. As expected p{Aa,T) > Poa{T) and p{Aa^T) > 
Poo{T) (thin lines in Fig. i.e. the ratio in question is 
> 1. Accounting for the complement via equations © 
and © (using values in Table P) collapses results from 
all the calculations to a single line recovering the bulk 
behavior (thick lines in Fig.^J; i.e. the ratio in question 
is 1. Deviations at low T are due to the increasing 
effects of monomers which have cq = 8 (12) in c? = 2 (3). 



In conclusion, a general approach in terms of the com- 
plement has been developed which allows one to account 
for finite liquid drop size effects and to extrapolate from 
finite to infinite systems. We have demonstrated the ap- 
plicability of our method using lattice gas model (Ising) 
calculations. This method can be generalized to include 
other energy factors present in the nuclear case (e.g. sym- 
metry. Coulomb and angular momentum energies) which 
is another important step towards determining the liquid- 
vapor phase boundary of infinite, symmetric nuclear mat- 
ter from experimental nuclear data. 
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